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Abstract. The conditions for emergence of Kolmogorov turbulence, and related weak wave turbulence, in 
finite size systems are analyzed by analytical methods and numerical simulations of simple models. The 
analogy between Kolmogorov energy flow from large to small spacial scales and conductivity in disordered 
solid state systems is proposed. It is argued that the Anderson localization can stop such an energy flow. 
The effects of nonlinear wave interactions on such a localization are analyzed. The results obtained for finite 
size system models show the existence of an effective chaos border between the Kolmogorov- Arnold-Moser 
(KAM) integrability at weak nonlinearity, when energy does not flow to small scales, and developed chaos 
regime emerging above this border with the Kolmogorov turbulent energy flow from large to small scales. 
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1 Introduction 



The concept of Kolmogorov turbulence 1941 [1,2] describes 
how the energy flows from large to small space scales in a 
turbulent regime. According to this concept an energy is 
injected on large scales, e.g. by wind, and it is absorbed 
on small scales due to dissipation. As a result a stationary 
algebraic distribution eu oc fc~ 5 / 3 of energy flow is estab- 
lished over wave modes k [1,2]. This concept was shown 
to be generic not only for hydrodynamic turbulence but 
also for other types of nonlinear waves. This phenomenon 
became known as the weak turbulence. With the help of 
diagrammatic technique Zakharov and Filonenko 1967 [3] 
derived the kinetic equation for weak turbulence of capil- 
lary waves and demonstrated the existence of stationary 
algebraic energy flows similar to those of Kolmogorov [1]. 
Later the concept of weak wave turbulence was generalized 
for various types of nonlinear waves as it is described in 
detail in [4,5]. However, an enigma of turbulence still re- 
mains as it is clearly stated in a poetic claim [6] "Through 
mechanisms still only partially understood, wind transfers 
energy and momentum to surface water waves. " 

Indeed, the kinetic equation for energy flow from small 
to large wave vectors is derived in the regime of weak 
nonlinearity and random phase approximation, known in 
plasma as a quasilinear approximation [7] . This hypothesis 
is at the basis of the whole theory as it is directly stated by 
Zhakharov and Filonenko at the first paragraph of their 
fundamental paper [3] : "In the theory of weak turbulence 
nonlinearity of waves is assumed to be small; this enables 
us, using the hypothesis of the random nature of the phase 
of individual waves, to obtain the kinetic equation for the 
mean square of the wave amplitudes". However, in finite 



size systems one has a a discrete spectrum of linear modes 
and the dynamical origins and conditions for validity of 
the random phase approximation and hence, for validity 
of the Kolmogorov concept of energy flow from small to 
large space scales, are still to be established. 

In fact, it is known that a weak nonlinearity does not 
always lead to global ergodicity and chaos over all wave 
vectors or over all degrees of freedom. The well known ex- 
ample is the Fermi-Pasta-Ulam (FPU) problem [8] where 
nonlinearity should be strong enough to generate global 
chaos as it was pointed by Chirikov and Izrailev [9]. At 
present the FPU problem still remains under active stud- 
ies aiming to understand its ergodic properties in the limit 
of low energy or weak nonlinearity [10,11]. For the FPU 
problem it was shown that, in the limit of very small 
nonlinearity combined with the resonant approximation, 
which is typically used for the derivation of the kinetic 
equation for wave [3-5], it is possible to have chaos only 
at small k— vectors with exponential decay of energies at 
high /c— vectors (see [12] and Fig. 2 therein). Thus this re- 
sult [12] shows that certain conditions are required for the 
emergence of algebraic stationary flows in the weak tur- 
bulence in finite systems. 

Indeed, in the case of a few degrees of freedom it is 
established that in the limit of weak nonlinear perturba- 
tion almost all phase space of a typical Hamiltonian non- 
linear system remains in an integrable regime known as 
the Kolmogorov-Arnold-Moser (KAM) integrability (see 
e.g. [13,14] and Refs. therein). The rigorous form of this 
statement is known as the KAM theory. The transition 
to global chaos and ergodicity requires that nonlinearity 
exceeds a chaos border which can be determined by nu- 
merical simulations or, in a number of cases, analytically 
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by the Chirikov criterion [13, 14], [15]. Below the chaos bor- 
der a main part of the phase space remains intcgrable and 
a chaotic spreading is possible only due to the Arnold dif- 
fusion via tiny chaotic web in a separatrix vicinity [13, 
14], [16]. In the limit of small nonlinear perturbation the 
measure of these chaotic layers drops exponentially [13, 
14] even if for a larger number of degrees of freedom this 
exponential regime can become valid only at very very 
small perturbations [16,17]. 

Of course, the initial concept of Kolmogorov turbu- 
lence and weak wave turbulence considers systems of large 
scale with a continuous spectrum of linear waves [4,5]. 
However, the modern experiments on wave turbulence are 
done with systems of finite size (see e.g. [18-20]). Also all 
numerical simulations are done with the finite size sys- 
tems (see e.g. [21, 22], [23], [24,25], [26,27]). It is recog- 
nized that in such finite systems certain resonant condi- 
tions play important role for energy flows in k— space [21, 
23], [24,25], [26]. Indeed, it is natural to assume that a 
discreetness of linear frequencies in a finite system can 
stop chaotic spreading if nonlinear frequency shifts be- 
come smaller than a typical spacing between linear reso- 
nant modes. Such a criterion has been put forward in [12] 
where its validity was confirmed for the FPU model in the 
resonant approximation. 

However, the situation may be even more complicated 
and the energy spreading can be suppressed even in sys- 
tems with everywhere dense spectrum. Indeed, in disor- 
dered systems it is known that the spreading of quantum 
(linear) waves can be stopped due to the phenomenon of 
Anderson localization [28,29] even if the spectrum of lin- 
ear waves is everywhere dense and classical particles can 
diffusively spread over the whole system. The effects of 
weak nonlinearity on the Anderson localization are now 
under active investigations of different groups (see e.g. [30, 
31], [32,33], [34,35], [36,37], [38,39], [40,41]). Usually the 
numerical simulations are done for a discrete Anderson 
nonlinear Shroringer equation (DANSE), which allows to 
perform numerical simulations in an efficient way up to 
very large times and thus to study wave packet spreading 
in space (see e.g. [33]). It is established that at moderate 
nonlinearity f3 a subdiffusive wave packet spreading con- 
tinues up to enormously large times being at least by 8 
orders of magnitude larger than a typical time scale in a 
system [30,33], [36,37], [39]. However, at small nonlinear- 
ity being below a certain border f3 c (f3 < j3 c ) the spread- 
ing is absent up to maximal numerically available times. It 
should be pointed out that the exact mathematical results 
are difficult to obtain even in the limit of f3 — > since the 
spectrum of the linear problem is everywhere dense so that 
resonances appear on large space scales [35] . Due to that 
it is difficult to develop the mathematical KAM theory in 
such a regime. 

Below the critical value of nonlinearity /? < f3 c the wave 
amplitudes decay exponentially inside a disordered layer 
[31,33], [34,40] in a way similar to a disordered linear me- 
dia in the regime of the Anderson insulator characterized 
by an absence of conductivity in such a system [29]. 



It is rather appealing to expect that the spreading in 
k— space of Kolmogorov turbulence will go in a way simi- 
lar to the DANSE case. In fact, for the quantum Chirikov 
standard map, known also as the kicker rotator, it is es- 
tablished that for a periodic driving in time the quantum 
localization of dynamical chaos takes place in the momen- 
tum k— space in a way similar to the Anderson localiza- 
tion in a coordinate space [42,43], [44]. The extension of 
this linear wave model to the case of the kicked nonlin- 
ear Schrodinger equation (KINSE) was proposed in [45]. 
In this KINSE model the nonlinear wave interaction takes 
place locally in space while we are interested in the en- 
ergy spreading in the momentum k— space as it is usually 
the case for the weak wave turbulence [4] . In this respect 
the situation is different compared to the DANSE model 
where both nonlinear wave interaction and spreading take 
place in coordinate space [33] . Even if it was argued that 
there is a certain similarity between these two cases [30] 
a special more detailed analysis of the KINSE model is 
required. In this work the KINSE model is studied on a 
large time scales and the links with the DANSE model 
are traced in a firmer way. The implications of the ob- 
tained results for the Kolmogorov turbulence in finite size 
systems are discussed. 



2 An example from the FPU problem 

Let us discuss briefly the effects of discreteness of lin- 
ear wave spectrum on example of the a— FPU problem 
following the results presented in [12]. It is shown there 
that in the long wave limit the system dynamics can be 
described by an effective renormalizcd Hamiltonian Hrm 
(see Eq.(4) in [12]). This Hamiltonian is only a resonance 
approximate description of the initial FPU problem. How- 
ever, it is important to see what are the properties of this 
Hamiltonian itself, since it has a typical form of resonant 
Hamiltonians considered in the theory of weak turbulence. 
This Hamiltonian Hrn = -Hflivo + HjiNi nt has an unper- 
turbed part Hrno oc k 3 Jk corresponding to the renor- 
malized linear spectrum of long waves and a part describ- 
ing the renormalized resonant interacting waves HRNint oc 
cos((f>k2+k 1 — 4>k 2 ~~ i'ki ) usually used in the theory of weak 
turbulence [4] (here (Jk,<pk are conjugated pairs of action- 
phase variables). It is shown in [12] that even if the dynam- 
ics of Hamiltonian Hun is chaotic for waves with k ~ 1 
this chaos does not create energy flow to high wave vectors 
k and the energy density drops exponentially at large k 
(see Fig. 2 in [12]). This shows that the random phase ap- 
proximation assumed in the weak turbulence theory can 
be not correct and that there can be no Kolmogorov flow 
from small to large wave vectors in finite size systems. 



3 The KINSE model description 

We focus here on the KINSE model described by 



iHd^/dr = -d 2 ^/2d 2 x + /3|</>| 2 V> 

-k COS X if) Em=-oo S ( T - mT ) ' 



(1) 



D.L.Shepelyansky: Kolmogorov turbulence, Anderson localization and KAM integrability 



where j3 characterizes nonlinearity, k is kick potential am- 
plitude, T is a period between 8— function kicks periodic 
in time, in the following we put H = 1. We impose the peri- 
odic boundary conditions with ip(x + 2Tr) = ip(x) and nor- 
malization condition L \ip(x)\ dx — 1. The linear wave 
expansion has the form ip(x) = ^2 n -0„exp(— inx)/\/2/K 
with conserved normalization ~^2 n \ip n \ 2 — 1- The second 
moment is defined as <j(t) — ^2(n — no) 2 \ip n (t)\ 2 , in the 
following we measure time t in the number of kicks (no is 
the initial mode). 

The linear model at j3 = has been studied in a great 
detail (see e.g. [42, 43], [44]). The semiclassical regime of 
this model corresponds to k ^> 1, T <C 1 with the classical 
chaos parameter K = kT = const. In fact K is the chaos 
parameter in the Chirikov standard map which describes 
the classical dynamics [13,14]. The classical dynamics is 
globally chaotic for K > 1 with a diffusive growth on en- 
ergy proportional to a with a = Dt and the diffusion rate 
D w fc 2 /2 for K > 4. In the quantum case this diffusion is 
localized due to quantum interference effects with expo- 
nentially localized Floquet states \ip n \ ex exp(— 2\n— no\/£) 
and the localization length £ « D/2 [42,44]. This dynam- 
ical, or Chirikov localization, is similar to the Anderson 
localization in disordered linear lattices where the mo- 
mentum states n play the role of spatial coordinate [43, 
44]. In a difference from the Anderson localization, which 
takes place in presence of disorder, the Chirikov localiza- 
tion takes place in a purely dynamical system without any 
randomness, but due to dynamical chaos the mechanism 
of localization of dynamical diffusion is similar to the one 
of Anderson localization. In absence of kick (K = 0) the 
model is reduced to the integrable nonlinear Schodinger 
equation. The KINSE model has been realized experimen- 
tally with the cold atoms and Bose-Einstein condensates 
(BEC) in kicked optical lattices (see e.g. [46,47]). The 
Chirikov localization at /3 = was observed experimen- 
tally [46], the studies of effects of BEC nonlinearity on 
this localization are now within experimental reach [47] . 

The KINSE model was introduced and studied in [45] . 
It was shown there that a narrow soliton has a long live 
time during which it follows an integrable or chaotic tra- 
jectory of the Chirikov standard map. In a regime when a 
soliton is destroyed it was found that there is still a sup- 
pression of classical diffusive growth of a. In [30] it was 
conjectured that this growth is similar to the case of the 
nonlinear kicked rotator (KNR) model, where there is a 
nonlinear phase shift for linear modes in momentum rep- 
resentation (it takes place during each kick period and is 
proportional to f3\if> n \ 2 ). The evolution of the KNR model 
is described by a nonlinear map for the wave function [30] : 

V>„ = exp(—iTn 2 /2 — i/3\ip n \ 2 ) exp(— ik cos x)ip n . (2) 

Here the bar marks the wave function after one period 
of perturbation, the operator cosf gives the Bessel cou- 
pling between the momentum states ip n . It was also argued 
[30] that the KNR behavior is similar to the case of the 
DANSE model, where the linear modes are also exponen- 
tially localized. The DANSE has the form 

idip n /dT = E n ip n + /3| ip n \ 2 ^ n + V{ip n+1 + tp n -i) , (3) 



where /3 characterizes nonlinearity, V = 1 is a hopping 
matrix element, on-site energies are randomly distributed 
in the range — W/2 < E n < W/2. The spreading in this 
model at moderate /? ~ 1 is characterized by a subdiffusivc 
growth 

a(t) ex t a (4) 

with the exponent a ss 0.3—0.4 (sec details in [33, 36], [37]). 
The similar values of the spreading exponent a have been 
found for the KNR model [30, 36], [48]. 

If to consider the lattice sites n in (3) as the momen- 
tum states then it is clear that the nonlinear coupling in 
(3) is local while in the KINSE model (1) it is strongly 
nonlocal in momentum or k— space, as it is usual for non- 
linear wave interaction in the regime of weak turbulence. 
Hence, the KINSE model is more adapted to the studies 
of the Kolmogorov turbulence in finite systems. Indeed, 
the kicks take place on a spatial scale of the whole sys- 
tem (of size 27r) and can be considered as a model of wind 
which pumps energy from small k— vectors (small |n|) to 
large ones. The numerical studies of the model (1) are 
presented in the next Section. 



4 Numerical results for KINSE model 

To integrate numerically the evolution described by (1) it 
is convenient to use unitary small step integrator making 
small kicks in coordinate space with the local space non- 
linearity and returning back and forth to the momentum 
space n with the fast Fourier transform as described in 
[45]. However, on very large times t ~ 10 6 the nonlinear- 
ity can generate exponential instability on high fc— modes 
which is of a purely numerical origin related to discretiza- 
tion. To eliminate this artificial instability the simulations 
are done on each step in an enlarged space size (approxi- 
mately 4 times larger than the size of physical modes) and 
after each small step the amplitudes out of the physical 
size are suppressed to zero. Such a method is similar to the 
aliasing approach [49] which efficiently suppresses numeri- 
cal instability on high modes [50] . Usually the simulations 
are done with the total number of states N to t — 2 10 and 
the size of physical state N = 220 (-110 < n < 110). 
The number of small steps N s on one kick period varied 
between 100 and 1000 with a special check that it does 
not affect the accuracy of the results. The numerical inte- 
gration method preserves the total probability up to a nu- 
merical double precision. The maximal value of t reached 
in the numerical simulations is t = 10 7 . 

The time evolution of the second moment o~(t) and 
probability distribution over linear modes w n = \ipn\ 2 in 
(1) are shown in Fig. 1 for the case of moderate nonlinear- 
ity /3 = 1 and a small kick amplitude k = 0.3 correspond- 
ing to the chaos parameter K = 0.6 being below the global 
chaos border K c w 1 for the Chirikov standard map [13]. 
It is clear that there are only quasi-periodic oscillation of 
a and that the distribution in momentum space remains 
exponentially localized. Thus, in this regime there is no 
energy flow to small scales and random phase approxima- 
tion assumed in weak turbulence [3,4], [5] is not valid. This 
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Fig. 1. Left panel: dependence of the second moment a on time 
t measured in number of kicks in (1). Right panel: probability 
distribution w n over linear wave modes n at times t = 10 3 
(blue), 10 4 (green), 10 5 (magenta), 10 6 (red) (all curves are 
superimposed). Here f} = 1, T = 2, k = 0.3, K = kT = 0.6 the 
initial state is at zero mode n = 0. 



i [44]). A very weak nonlinearity j3 = 0.05 does not af- 
fect this localization which persists up to maximal times 
t = 10 7 reached in numerical simulations. We note that 
following [33] the values of a are averaged over logarithmi- 
cally equidistant time intervals that suppress fluctuations 
at large times. 

However, at moderate nonlinearity (3 = 0.5 and (3 = 1 
the second moment shows a subdiffusive growth with the 
algebraic exponent a ss 0.4. The exact fit values of a are 
given in the caption of Fig. 2. As in [33] the statistical error 
bars are relatively small but there are long time correla- 
tions that can affect the real a value on large time scales 
as discussed in [33,36], [37,38]. A similar value of a « 0.4 
has been found by A.S.Pikovsly for numerical simulations 
of the KINSE model by another numerical method for 
t < 10 6 [51]. 



result is similar to a usual observation that a small wind 
(small k here) is not able to produce a turbulent storm. 
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Fig. 2. Dependence of the second moment a of probability 
distribution over wave modes n on time t for (1). Here, k — 3, 
T = 2, K = kT = 6, P = 0.5 (blue curve) and /3 = 1 (red 
curve); a(t) is averaged over logarithmically equidistant time 
intervals; the initial state is no = 0. The dashed line shows an 
anomalous diffusion a ~ t a with the exponent a = 0.4; the fit 
of data in the range 3.5 < log 10 1 < 7 gives a — 0.346 ± 0.014 
(for /3 = 0.5), a = 0.438 ± 0.007 (for j3 = 1). The horizontal 
dashed black curve shows data at j3 — 0, it practically coincides 
with the data for j3 = 0.05 shown by green curve. 



Let us now consider the regime above the classical 
chaos border with K = 6. In his regime the classical sys- 
tem, described by the Chirikov standard map [13], has 
diffusive energy growth with a sw k 2 t/2 w At. The re- 
sults of Fig. 2 at (3 — show that this diffusion is local- 
ized by quantum interference with the localization length 
£ rj fc 2 /4 w 2.2 and er ~ £ 2 ~ 10 (the value of sigma is 
slightly higher than £ 2 due to mesoscopic fluctuations of 
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Fig. 3. Probability distribution w n over linear wave modes n 
in (1) at times t = 10 3 (blue), 10 5 (green), 10 7 (red) for = 0.5 
(left panel), 1 (right panel); other parameters are as in Fig. 2. 
The dashed curve shows the probability distribution at /? = 0, 
t= 10 7 . 




Fig. 4. Probability distribution w n over linear wave modes n 
in (1) at time 10 7 for /3 = 0.0 (dashed black), 0.05 (blue), 0.5 
(magenta), 1 (red); other parameters are as in Fig. 2. 
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The evolution of probability distribution w n with time 
is shown in Fig. 3 for j3 — 0.5; 1. It is clear that nonlin- 
earity destroys localization: an approximately flat plateau 
of probability is formed ( "chapeau" ) which size is slowly 
growing with time. The data for w n at large time t = 10 7 
show that there is a significant increase of the distribution 
size for j3 = 0.5, 1, while for /3 = 0.05 the probability re- 
mains localized in a way similar to the linear case [3 — 0. 
These data qualitatively confirm existence of chaos bor- 
der in nonlinearity with a certain j3 c ~ 1/10. However, it 
is not excluded that some very slow processes related to 
the Arnold diffusion can lead to "escape" of some small 
probability to larger n values at exponentially large times. 

The results obtained for the KINSE model (1) show 
that the behavior of this model is qualitatively similar 
to the one found in the KNR (2) and DANSE (3) mod- 
els: the spreading in the momentum remains localized be- 
low a certain chaos border in nonlinearity strength, while 
above this border the spreading continues in a subdiffu- 
sive way with the algebraic exponent being close to the 
value a ps 0.4 found previously for KNR and DANSE. The 
main difference from DANSE is that for KINSE we have 
interchange between momentum and coordinate spaces. 
However, such a similarity between momentum and co- 
ordinate space is well known [43,44] and was already dis- 
cussed for the KNR model in [30] . The main new aspect of 
the KINSE model is a long range coupling in the momen- 
tum space due to the locality of wave interaction in the 
coordinate space. However, due to a localized nature of 
linear eigenmodes in the momentum space this long range 
coupling still can give transitions only on a size of local- 
ization length I and this does not produce a qualitative 
difference from the case of short range interation appear- 
ing in DANSE. A similar situation appears for the model 
of two particles with Coulomb interaction in the regime of 
Anderson localization as it is discussed in [52]. 



5 Kolmogorov turbulence in Sinai billiard 




Fig. 5. Dependence of the second moment a in the DANSE 
model (3) at W = 4, V = 1, /3 = 1 (red curve) and in the 
DANSE model with additional static field potential (SDANS- 
BIL model)) with SE n = f\n\ at / = 0.5 (blue curve) at the 
above parameters and the same disorder realization. The initial 
state is at n = 0. 





-1000 -500 500 1000 



Fig. 6. Probability distribution w n over linear wave modes n at 
times 10 8 for the two cases of Fig. 5 with the same attribution 
of colors; two panels show the same data on different scale. 



On a first glace one can get an impression that the KINSE 
model is a rather specific one: it is one-dimensional, there 
are kicks etc. However, it is known that the Chirikov stan- 
dard map is generic and describes a variety of real physical 
systems [13,14]. Also it is known that the Chirikov local- 
ization found first in the linear KINSE and NKR models 
aX (3 — Q appears in a variety of real systems (see e.g. [53, 
54]). Thus we expect that the behavior found here for the 
KINSE model at finite nonlcncarity j3 will appear in more 
realistic systems. 

Let us discuss an example of such a system. For that 
we consider the nonlinear Schodinger equation in a two- 
dimensional chaotic billiard (e.g. the Sinai billiard [55]) in 
a presence of a monochromatic driving: 

idip/dr = -Aip/2 + V(x, y)ip + /3\ip\ 2 ^ + F sin(uT)xip . 

Here the potential V(x, y) is determined by a rigid bound- 
ary of the billiard, F and u> are the amplitude and fre- 
quency of a monochromatic force acting on a particle in- 



side the chaotic billiard. We assume that the classical dy- 
namics inside the billiard is chaotic that leads to quantum 
chaos and ergodicity of the eigenstates of the linear prob- 
lem (see e.g. [56]). The microwave monochromatic force 
acting on a classical particle creates a diffusive growth of 
its energy with time but the quantum interference effects 
lead to exponential localization of this diffusion (at j3 = 0) 
[57], in a way similar to the KINSE and KNR models at 
/? — 0. It is interesting to know how the nonlinear term 
with (3 in (5) will affect this localization and if nonlinear- 
ity can create the Kolmogorov flow of energy from small 
to large wave vectors. An analogy with the models of the 
previous Sections gives an idea that at weak F and /3 there 
will be no energy flow to high level numbers of the linear 
quantum billiard. However, for the quantum billiard there 
is a new element which we discuss below. 

Indeed, the density of energy levels inside the 2D quan- 
tum billiard is approximately constant, up to quantum 
fluctuations, [56]. Due to that the energies of quantum 
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levels grow with the level number n, approximately as 
E n ~ pn where p is the average density of levels. Such 
a behavior of energy levels is different from the case of 
DANSE where all unperturbed energy levels are located 
inside a finite energy band. A linear growth of E n with n 
corresponds to a presence of a static Stark field with an 
additional term 5E n — fn. Such a DANSE model with 
a Stark field has been studied in [58] and it was shown 
that a subdiffusive spreading goes in a way similar to the 
DANSE model at moderate values of /. However, for a bil- 
liard there is a minimal energy so that we have only n > 
(like a triangular potential). To model such a situation one 
needs to assume that in the DANSE model (3) there is 
an additional energy shift SE n = f\n\. Thus the DANSE 
model with such a modulus Stark term reproduces the 
energy growth with level number typical of the quantum 
billiard. We will call this the Stark DANSE model of bil- 
liard (SDANSEBIL model). The numerical simulations of 
this SDANSEBIL can be done in the same rather efficient 
way as it is described in [58]. The comparison of the re- 
sults at moderate nonlinearity /3 = 1 for / = (DANSE) 
and / = 0.5 (SDANSEBIL) is presented in Fig. 5 and 
Fig. 6. It shows that a finite / reduces the value of a al- 
most by two orders of magnitude and that the probability 
distribution over levels is localized in a much stronger way 
compared to the case with / = 0. The physical origin of 
this suppression should be attributed to the energy con- 
servation which for / > leads to a more rapid decrease of 
probability on high levels: for a homogeneous probability 
distribution on levels < n < n max the energy conser- 
vation imposes \if) nmalc \ 2 ~ l/{f n max) while at / = the 
restriction from norm conservation gives slower decay of 
probability with |VVi ma J 2 ~ l/n max . Due to this reason 
the spearing over levels is suppressed stronger in the pres- 
ence of finite static field / and triangular form of energies 
5E n oc \n\. We note that the case with SE n ex n, con- 
sidered in [58], does not give additional restrictions due 
to cancellations of negative and positive n contributions. 
The obtained results show that the transition to ergod- 
icity in the SDANSEBIL model is practically absent and 
thus there is no energy flow to high wave vectors. 

Due to similarity between DANSE and KINSE and 
KNR models discussed in previous Sections we make a 
conjecture that the behavior similar to the one found for 
SDANSEBI1 model will take place for the evolution in the 
nonlinear billiard model (5). As a results of that observa- 
tion it is possible to make a conjecture that there will be 
no energy flow to high wave vectors for the Kolmogorov 
turbulence in Sinai billiard described by Eq.(5). Of course, 
it would be very interesting to perform direct numerical 
simulations of the model (5) but this would require much 
more advanced and heavy numerical simulations. 



6 Discussion 

In this work we discussed the properties of weak wave 
turbulence in finite systems. On the basis of numerical 
simulations and analytical results we argue that the dis- 
crete spectrum of linear frequencies, typical for finite sys- 



tems, imposes specific conditions for appearance of the 
Kolmogorov energy flow from large to small spacial scales. 
In absence of random phase approximation such a flow can 
appear only above a chaos border at a sufficiently large 
nonlinear coupling between linear modes and/or strong 
driving force. The considered models show that the An- 
derson localization, appearing in the wave vector space of 
linear modes, can stop the energy flow from large to small 
scales if nonlinearity is below the chaos border and the sys- 
tem is in the regime of KAM integrability. A similar situa- 
tion appears at a small amplitude of energy pumping. In a 
qualitative way such a regime corresponds to a small wind 
which cannot generate turbulent ocean waves. Of course, 
the models analyzed here are relatively simple and hence, 
the numerical and experimental studies of more realistic 
systems, like e.g. the model (5), are highly desirable. Such 
studies will allow to understand a new regime of nonlin- 
ear waves where interplay between Anderson localization, 
nonlinearity KAM integrability and Kolmogorov turbu- 
lence in finite systems opens new interesting and unsolved 
questions. 

Finally, it is useful to note that a quasiperiodic driving 
of the KINSE model at j3 = with two incommensurate 
frequencies can create the Anderson transition with ap- 
pearance of energy flow to high momentum states [59]. 
Such a transition has been observed recently in experi- 
ments with cold atoms in kicked optical lattices [60] . Thus 
it is possible that a transition to weak turbulence in finite 
systems is somewhat similar to the Anderson transition in 
disordered solids. In such a scenario the regime of Ander- 
son insulator corresponds to laminary waves and absence 
of energy flow from small wave vectors to large ones, while 
the metallic phase allows to have a turbulent energy flow 
from small to larger wave vectors. The analysis of possible 
links between these phenomena requires further investi- 
gations. The modern techniques of ultra cold atoms and 
BEC (see e.g. [46, 47], [60]) allow to study experimentally 
the effects of nonlinear wave interactions in such systems. 

I thank A.S.Pikovsky for stimulating discussions and 
critical remarks. 
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